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ABSTRACT 


We consider the optimum design of ph° to ^ ^ mht condit i 0 ns. Extremely compact Q 
airborne and spacebome platforms under berth y JP ^ ^ rates and can 

switched microlaser transmitters P rodu ® Drec ise ranging. To guide the design, we have 

easily generate subnanosecond pu sew simple algorithms, based on Post-Detection 

modeled the solar noise background and devel< op P J , from a high no ise 

Poisson Filtering (PDPF), to optimally as detector and/or 

background during daytime operations. analysis We describe an airborne 

receiver dead times, have also C0 “^^nent Incubator Program, which is designed to 
prototype, being developed under NASA s Instrum energies on 

Operate at a 10 kHz rate from aircraft cruise aft efficlgnt s F ystem designe d to 

order of a few microjoules. W ® a ‘ s0 Sample the Martian surface at rates up to 4 .3 

°^r, 0 SSwXtr^t' g a l iS roughly 400 dmes greater (4 kHz vs 10 Hz). 
Relative to conventional high power laser 

altimeters include: (1) a more efficien u laser powe r-telescope aperture product, 

of magnitude greater surface sairjdmg ; rates for ^ ^ ^ vo|um< _ and weight of the 
(2) a simultaneous two order of resolve the source of the surface return in a 

telescope system; (3) the unique ‘ability to or imaging detectors; and (4) improved 

photon counting mode through the us P from feoth (1) md (3 ). Furthermore, because o 

vertical and transverse spatial resolutio • ® himeter is inherently more eyesafe to 

r~~d »dLrpSnete mten,a, optical damage, which can tennhiate a space 

mission prematurely. 
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1. INTRODUCTION 

Spacebome laser altimeters typically use modest energy (50 to 100 mJ) . soli f 
telescopes (50 to 100 cm diameter), and high detection thresholds to achieve unambitious 
surface returns with few or no “false alarms” resulting from solar background noise [Abshire et 
al 2000- Dubayah et al, 1997; Ramos-Izquierdo et al, 1994], As a result of this convention 
design philosophy spacecraft prime power and weight constraints typically restnrt operations 

the order of a few tens of Hz which at typical planetary <** 9 ^ 
velocities of a few Km/sec, limits along-track spatial samplmg to one sample every 
meters There is great scientific interest in obtaining higher along-track resolution and/or bett 

cr^krack coverage, bu, ach,e™g .his capability toougb - 7 ". ^ 

(power) is not practical. This is especially true of altimeters destined for orbit about distant 
planets where instrument mass and prime power usage is severely constrained. 

We will demonstrate in this paper that the conventional high surfece 

anoroach to laser altimetry does not make efficient use of the available laser photons. The surface 

raun. rate of an oAkin^ritimeter can be increased by up to two 

laser output power by emitting the available photons m a high frequency (few KHz) tram ot low 
energy (< 1 mJ) pulses, as opposed to a low frequency tram of high energy pulses, and y 
employing single photon detection [Degnan and McGarry, 1998]. This mode °P era *° n 
2'Tdte risk of internal oprical damage tothe laser, Aereby improvntg 

and makes the beam inherently more eyesafeto a ground-based observer. In addrtion, Aese high 
be accomplislied wkh much smaller telescope apertures on Ae order of 10 to 20 
cm diameter Indeed, the contrast of the terrain “signal” against the solar-induced noise 
background is actually enhanced through the use of small receive telescopes As we will 
demonstrate, relatively simple onboard software algorithms, loose ^ 

Poisson filtering techniques, previously used in lunar laser t*bb<* Ze toiafa J ’ 

can be employed to identify and extract the surface samplmg data from the more abuadan 

noise prior to onboard data storage or transmission of the data to a ground station. 

The near order of magnitude reduction in telescope diameter greatly simplifies the mechanics of 

scanned systems and allows the use of relatively inexpensive, modest diameter optica we ges 

holographic optical elements to simultaneously scan both the 

cross -track interrogation of the terrain while still maintaining ^ 

background noise suppression. In addition, since the volume and weight of a telescope an 

support stturture noSnally varies as D 2 where D is the telescope diameter a 

twcforder of magnitude reduction in weight and volume and a comparable reductio 

fabrication cost can be realized. 

Further unique performance enhancements are also possible when operating in a P hoton ^ 0 ^tmg 
mcT ZZercial and developmental photon-counting detectors now exist wh.ch are capable of 

ranging resoLon as well as 

photon event within the receiver field-of-viev. These devices * 

designed high bandwidth photomultipliers [Charboneau et a , , , ’ . Q q 71 

Sorsky^t al 1996] or as Avalanche PhotoDiode (APD) arrays [Vasile et al, 1997], Wdh 
high angular resolution of the single photon source, the measured range becomes near y p 
3omt measurement, i.e. from the internal altimeter reference point to a small an* of 
uncertainty on the surface, which is determined by the angular resolving power of the 
phcZSor and can be very small compared to the total beam area on the ground. The ranging 
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precision is now limited by the laser pnlsewidUr, the hmmg capab.hnes of Ure , range 
tiie much-reduced residual spreading of the return waveform , or Poisson generating 

rr” P h^ 

ss 

We believe that the S 

al, 1997] It is also an attract, ve successor u mappmg the surface 

fMOI instrument rRamos-Izquierdo et al, 1994J, wnicn is y yy /\ArQ\ 

instrument sensitivity is largely Jlftabase of land surfaces 

n to 12 km) Even with their higher pulse energies, conventional airborne laser attimet 
^pT ca Uy must X* altitudes below 3 km, which usually requires a special waiver from 
government aircraft authorities, 

hr Section 2 and the appendices, we develop the *f°? 

Section 3 provides a brief technical description an p Incubator Program. Section 4 

provides a brief summary of the results. 

2. THEORY OF OPERATION 

2.1 PHYSICAL MOTIVATION FOR THE "MICROALTIMETER" 

altimeter equation [Stitch, 1972], 
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r r 2l^M _ r M. =r P & rr ^ 

n s =[pcosaT Q \-^r-t a r2 a 

, r • t . p altimeter constant h v is the laser photon energy, r\ q is the detector quantum 
OT^wayattnospheric transmission along the altimeter nad,r pad, graen by 


T 0 = exp 


( 00 ^ 
-jdzju(z) 

V * 


(lb) 


where Ms the height of the surface at»ve sw, level ******* 
atmospheric extinction coefficient. Solvmg for fgs m ( ) yi 


V q Vr P T A R 
\ s nhvR 2 


( 2 ) 


f QS =[P cosaT *]~ n 

Poisson process, i.e. 


if 

/r - fQS^D{ n s’ n t) “ /max ^ ^ “ A 


r=o^’-y 


(3) 


where »,is the receiver detection threshold in photoelectrons, and the maximum surface return 
rate is defined by 


%T]A^r 

nhvR 2 


(4) 


f^.=[p cosdr «\ 

where we have divided die latter expression 

and a term which depends solely on instrumen nnerture nroduct the maximum 

Figure 1 demonstrates that, for a given laser power rtect j on threshold («r=l) and a low 

surface return raters obtamed by plies using a high repetition rate, but low 

mean signal strength, i.e. n s < <->• Hnysiw y, v ^ must be taken mto account in 

energy, laser transmitter, but cl ^ rl >^ e ® m X Editions from orbit. Our prototype 

systems designed to operate under both y use eated visible detectors, which are 

^Ty^^ 

of — noise dunng daylighi 

operations. 
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2.2 CORRELATION RANGE RECEIVER 

A correlation range receiver (CRR) can be need to distinguish weak sigwls ^ m ^ cb ^ h « on 
noise counts during daylight operations (Degnan and Mctony j W Jhe CRKmeas ^ 
time-of-flight at multi-kilohertz rates rdative to the outgoing p w m flight with 

high repetition only fcie photons which e£er the receiver dunng a 

detoe^by the 

occurring wither a frame as defined by 2 is called a cell A set 

of 10 s -“ 

The terrairfis shown as a solid black line passing through the various cells. 

Software mthe CRRmntatively . dentrfies 

STfs^eVS'on 2 3 5). tf the count '““^^'^itively [Sfi^SStoin^ 

7, and 8. It is also statistically possible that a noise ^ ^ confinition of signal acquisition 

signal. Under these weaker signal to noise co > rThterton et al 1998] to the data which 

can be obtained, if necessary, by applying an of |j (j) passes the threshold 

in 7 of 10 frames, it also illustramd in F.gure 

2 . 

More sophisticamdtiackingalgorfthrns, which 

are also possible and can further unprove ?C!ri in time and usmg our 

terrain features and recover missing data. Forexample, 8 in bin 4 5 or 6 (cells 

adjacent range bin criteria, the lost signal in ^ ^ igna j ^ frame 2. However, looking 

^pmgt ”f ribmt Slhatrdrwhrchmreahtycomams most ofthe signal. 

AS a second example, we cons.der Frames 7 rmd 8 whereftre 

consecutive frames. Looking forward in h me ° Lookimfbackwaids in time from Frame 9, 
through 8 in frame 7 and bms 5 7. The overlapping 

however, suggests bins 4 through 6 in Frame 8 and a tnrougn contain the 

(cross-hatched) bins are 6 and 7 in Frame 7 and bms 5 . mU m Frame 8^ ^ ^ 

missing signal counts. While very ab ~P' violent onTbin criteria, we will 
signal cell to jump many range bins in a single frame and vioiare o _ 
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demonstrate in later simulations that the algorithm quickly identifies the surface at the new 
elevation. 


Surfece slope and/or roughness within the interrogated surface area. A,, will produce a return 
waveform (or in the case of photon counting, a Poisson generating function) that is significantly 
broadened relative to the initial laser pulsewidth, which is typically on the order of a nanosecond 
DSributedCets, such as clouds or ground vegetation, also produce broadened r^ms and, with 
a properly designed receiver, are detectable usmg photon counting instalments [Ho and 
Pnedhorsky, 1995]. Multiple returns from different entities - e g. cloudtops , tree canopies 
and/or groiLd returns - can also occur on a single laser pulse, potentially causing several CRR 
cells, widely separated in range , to record valid returns. 


2.3 SOLAR NOISE BACKGROUND 

From equations (A7) and (A9) in Appendix A, the mean number of background noise counts 
generated within a given range bin by scattered solar radiation is given y 


n b =pA r z b =[Nl{M)Cl r T b _ 


VejjA 

7thv 


pT 0 Usece ‘ cosy/ + 


1 -T n 


1+sec 9 S 


4(l + sec0,) 


(5) 


where the first and second terms respectively correspond to the noise counts per range bin idue to 
scattered solar radiation from the surface under study and the intervening atmosphere. In Eq. (5) 
which is a generalization of Eq. (18) in [Stitch, 1972], P is the combined noise count rate per unit 
area of receive aperture, N° x is the wavelength-dependent solar spectral illuminance as measured 

just outside the planetary atmosphere (N x °=02 watts/m 2 -A° at 532 nm for Earth), AX ' js » the 
spectral filter bandpass, A is the the receiver field of view (FOV) m steradians^r, u Mto dur^on 
afthe range bin, yr is the subtended angle between the Sun and the surfece normal ^d ftisthe 
solar zenith angle at the time of the measurement. It is also shown in Appendix A tha 

* (6a) 

cosy/ — cos o’ cos 0 S + sin a sin 6^ cos^? 

where cp is the angle between the plane containing the altimeter line-of-sight and the Sim (x-z 
plane) Ld to ptae contoring to altimeter lineK,f-sigh. and to surface normal as m F.gnre 3. 

Since 1 > cos#> > -1 , we have the condition 


cos(o--^)>cos^>cos(<r + 0 S ) 


(6b) 


where cost has to maximum value of 1 when 9 = 0 and cr = ft and a mmimum value of 0 when 
the surface is in shadow (cos y/ < 0). 


2.4 INSTRUMENT CONTRAST AND SCALE ENERGY 
Using (1) and (5), we can define an instrument contrast 


W =1+ 4V!l,i^ = i4 

N b f Q s r f n b n » E sc 


(7a) 
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where we have assumed single photoelecton detection, N f -/«**/» the number of laser fires per 
frame and n and n h are the mean signal and noise counts per range bin per laser tire 
respectively. 'The contrast is equal to the ratio of total mean counts in the signal cell to the mean 
background counts collected within any other noise cell in the CRR We have also define a sc 

energy given by 


(7b) 




[A^”(AA)fi r /? 2r 6 sec ex] 

‘ cosv, +w |+sec ^ j 


as that transmitted single pulse energy which produces thesame mean signal count m a P£P^y 
sized cell (i.e. one containing all of the expected mean surface returns within a feme > 50 

noise background produces in any other arbitrary cell and which in turn producesa cell contrast 
of 2 The scale energy is independent of other details related to the instrument, atinosphere, or 
suic^properties siL these will have an equal scaling effect on both signal and background 
noise and hence will not affect the contrast. We note from (7b) that the scate ^rgyva 
quadratically with range and linearly with both the receiver solid angle FOV and the range bin 

duration. 

In the limit as T 0 ->1 (relevant to an altimeter orbiting an airless moon or asteroid), (7b) reduces 
to the simple expression 


E = N°,AAn r R 2 T h sec a cosy/ 1 < secrr 

SC A ' u 


( 8 ) 


Assuming a ground spot diameter of ~25 meters (corresponding to a 10 arcsecond laser beam at 
an orbitaf althude of 500 Km), a modest filter bandpass of ~3A°, a range bm of~200 ns^ an 
surface slope (o=0), the scale energy is less than 120 pJ, or almost 3 ordera of sn ^" er 

than the 100 mJ GLAS transmitter pulse energy. For a relatively steep surface slop , 

scale energy is increased by roughly a factor of 3 to 360 pJ. We will demonstrate shortly&ata 
transmittecTenergy 2 to 3 times the scale energy is usually adequate to a^eve high reha y 
discrimination of the surface returns from the worst case solar background. Thus, Earth mapping 
missions using less than 1 mJ of transmitted energy are feasible m photon counting mod 

In order to compute the maximum scale energy in the presence of a scattering atmosphere, we 
will find it useful to introduce the variables 


a) - 


1 

r»t + sec# r 


1 


< (jr) < CC 


(9a) 


and 


a = 


1 

4 p cos if/ 



which, when substituted into (7b), yield 


(9b) 
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E. = 


N°,A^ r R 2 r b sec a 


1 + a 


T 2 

i o 


COS^- 


ft> - 1 

ln<y 


( 10 ) 


cy 


Setting the derivative of (10) with respect to co equal to zero 

[-In 2 co + a(ln<y-<y + 1)1 

cos \ i/{ — rri 

r co In co 


dE N°,AAQ r R 2 T b sec<j 


= 0 


01 ) 


doo 


suggests a single extremum 


at ® =1. However, since E K and its derivative are both indeterminate. 


we 


must make repeated use of L’Hosprtal's rule to obtain 


lim E = 

m-*\ S 


and 


N°,AA£l r R 2 T b sec o~ 


£ 


1 

f i Y 

COS W -1 lit 

~Tp 

Ap 

IU 


(12a) 


dE N\A?El r R 2 t b sec a 


lim 
»-»> dco 


cos ^ ln li 


(12b) 


Since the derivative (12b) is negative the ^xtremum at a a ^ “ 

= 1 /T 0 2 , the smallest possible physical value of co in (9a), into (1U) yieias an pf 

scale energy, i.e. 


E < NlAAri r R 2 T b seccr 


1 

1 + — 
8 /? 


'±-1 

. 0 > 


(13) 


" Sr' 


low reflectivity surfaces that satisfy the condition 

532 nm (T 0 > 0.7), primarily affects surfaces with 

2,5 DISCRIMINATING S.GNALEROM NO.se WITH POST-DETECnON PO.SSON ANALYSIS 
To achieve a maximum surface return 

be so gr^Taf to! Ae p'uTse repetition rate (and hence surface sampling rate) for 


which, for typical atmospheric transmissions at : 
reflectivities less than about 0.13. 
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a given laser output power. In any case, we are interested in designing a CRR which can 
effectively identify signal cells under conditions of relatively low or modest contrast (C > 3). 


Since the system parameters must be chosen such that the mean noise count per range bin per 
laser fire is significantly less than the mean signal count per laser fire (which is also on the order 
of unity or less for efficiency reasons as previously described in Section 2.1), usually no more 
than one photon count (signal or noise) is accumulated in any given cell per laser fire. On the 
other hand, the total number of noise counts observed within the range gate per laser fire may be 
significantly greater than one since, although the probability of detecting a noise count in a given 
range bin is small, there can be a large number (hundreds to thousands) of range bins if the range 
gate (range uncertainty) is large. Furthermore, the requirement that the noise count per bin be 
much less than the mean signal count (<1) places a constraint on the maximum telescope aperture 
that can be employed in daylight operations to reach the necessary power-aperture product for a 
given surface sampling rate, although this limit can be relaxed somewhat by improved spatial and 
spectral filtering. 


From Poisson statistics, the probability that we will correctly identify and acquire the signal cell 
in a given frame is 


*-l M k 

P^ = P{N l >K) = l-e~ N ‘Y j - ‘ 


acq 


k~- 0 


k\ 


l-e 


CN, 


'MY 


k = 0 


k\ 


(15a) 


where N t = N, +N b = CN b is the total mean photoelectron count, equal to the sum of the mean 
signal and background counts in the cell, and K is the frame count threshold. Similarly, the 
probability of falsely identifying any given noise cell as “signal” is equal to 







(15b) 


Now, the mean number of noise cells per frame falsely registering a signal is given by 

^ false ^ birF false 


where N bm = r g / r b is the number of range bins within the range gate. Thus, for the M frames 
making up a superframe, we can define a differential cell count 


N(K) = M[P acq -N bm P false ] = M 


l-e 


-car, 


^(CN„)‘ 


k=0 


k\ 


-N, 


bin 


Jfc=0 * ! 


(17) 


which subtracts the expected mean number of falsely acquired noise cells from the expected mean 
number of true signal cell acquisitions occurring within a superframe. 

Choosing a frame threshold K in (17) that is too small will result in too many noise cells being 
falsely identified as signal and will drive the differential cell count to small («M), or even 
negative, values. On the other hand, choosing a threshold that is too high will reject many valid 
signal cells, along with the noise cells, and drive the differential cell count toward zero. We 
therefore pose the question - what integer frame count threshold, K, maximizes the differential 


9 



ceU count’ To ottarn lb. answer, we perform a quasi-differentiation whh respect to the frame 
threshold, K, on (17), i.e. 

(18) 


A K K 

Setting the latter expresston equal to aero and solving for K gives a simple formula for the 
optimum frame threshold , K 0 pt 

(C- lR+ln(/Q W,+ln(JV»,.) 

opt ~ In C InC 


(19) 


K, 


i.e. 


A 2 # 

AK 


wr n PO) 

2 = N(K opl + 1) - 2N(K opt ) + N(K opt - 1) - l C ) 


One can now use (.9) m (.7) to e, unmate *.* «. expmssion for the cptmuzed thfferentta, cel, 
count and obtain 


N op t =M\ 


1 -C «e 


(C- \)N„ _ 


-CN„ 


(c*-c^) 


m 


k=0 


k\ 


(21) 


Using (21), we can 


lim N a p = 0 

C-fl 


further show that, in the extreme contrast limits. 


(22a) 


and 


lim N o0 =M 

C-+* 


22b) 


Thus, since N opt is always positive for the and hence the 

identified signal cells always exceeds since that is the maximum 

latter can never exceed a mean ofone P er * theootimum threshold (19) into (15) and (16) 
number of allowed signal cells . One can P probability of falsely 

to compute the probability of corrertly F ftame which exceed the 

identifying a no.se cell as signa or the mean ^^ e Zt Zyses of airborne and 
optimized threshold. As we will demonstrate cell count arbitrarily close to the 

SS^emTre M £L* itfM si^ra, culls in a superfmme and zero false 
acquiskions) by choosing a sufficiently long frame mterval. 
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mean number of signal counts m the (*) '" ^ at higher con trast, fewer signal 

four values of the contrast, i.e. C 3, , . , ■> j re( j uce the number of false 

counts per frame are required to positively identify the si^al celled Into* a* ^ ^ 

alarms to an insignificant level. We fiirtl l^ differcr^Si cell Sunt and probability of correctly 
contrast values considered, the norma lz Irnm nnitv and the mean number of false 

per ^ Pro ^ e^, 

discrimination over a wide range of signal cell contrasts ( - >■ 

2.6 CHOOSING THE RANGE BIN AND FRAME INTERVAL 
The choice of range bin and frame interval are 

cannot be made arbtaarily large (e g., over a doped surface) wrthort evattua ly s e 

range bin . in order ,o rttain aU of the £££&£ the cell 

(10), a larger range bin also implies prey omonaUy larger sca^energje^ ^ ^ off , 

contrast for a fixed mean signal strength^ a ^° r frame ^ Q f the signal counts observed during 

constant steep slope o for the entire duraHon ’ chooseTrange bin duration which is 

the frame will always be contained within a single cell if we choose rang 

twice as large as the total signal spread within a frame 


r fc =2 


r»+-(D ! + ''?*r) ana 


= 2r, +Y r f 


(23a) 


where c is the speed of light, D, is the beam diameter on^egromrd, e 

according to ihe e.at.oo 


r s = r 0 + 


2D tan cr 


(23b) 


The parameter 


r = 


4 vf tan cr _ 4(v„ +2^/ Jcan tan^)tana 


(23c) 


is the rate at which the range bin grows with ■ ^hy“' Sret^mi^ual 

surface slope effects. For an unscanned system, *• ^ y for 

to the aircraft or ^ a ^ aft ° C ^’ ( | 3) suggestTthat the range bin duration can be 

1 forai1 butneariy vert,cai 


n 


slopes (<t « 90°). For typical aircraft cruise velocities on the order of 0. 1 to 0.2 km/sec, the range 
bin duration can be 5 to 6 orders of magnitude smaller than the frame duration. 

If we include a scanning mechanism on the altimeter to obtain cross-track measurements of the 
underlying terrain, the effective ground velocity of the beam can be significantly larger than the 
ground velocity of the host vehicle, leading to an even larger range bin and a further reduction in 
signal cell contrast. A relatively simple, although not necessarily optimum, scanning device is a 
slightly wedged optical plate spun at a frequency f scan in front of the transmit and receive optics. 
Because of the small optics (typically 10 to 20 cm) required for the microaltimeter, this is a viable 
scanning approach and produces a circular ground pattern superimposed on the linear motion of 
the aircraft or spacecraft. Since the frame interval is usually short compared to the circular scan 
period, the maximum effective ground velocity occurs when the linear ground velocity of the 
aircraft or spacecraft and the beam velocity due to the scan add vectorially as in (23c), where R is 
again the range to the surface and £ is the angular deviation of the circularly-scanned beam from 
nadir (half angle of the cone). Two counter-rotating optical wedges in series can be used to create 
a variety of scan patterns, including a linear raster scan in the cross-track direction (same rotation 
frequency and wedge angle but starting 180° out of phase in the along track direction). 


2.7 WAVELENGTH FIGURES OF MERIT 

Relative to 532 nm, the 1 064 nm fundamental Nd: YAG wavelength is known to have many 
"natural" advantages for altimetry, such as generally higher reflectance from various types of soil 
and vegetation, improved atmospheric transmission, and reduced solar background noise. 
Furthermore, one pays an additional penalty (~50%) in converting infrared photons to the visible 
via passive nonlinear crystals. However, it will now be shown, via two instrument figures of merit 
(FOMs), that the much higher detector quantum efficiencies and narrower, high throughput 
spectral filters available at 532 nm often makes the latter wavelength a better choice for the 
microaltimeter. 


X(nm) 

7* 

To 1 

7m 

P 

soil 

P 

dry veg 

P 

green 

veg 

FOMru 

soil 

FOMru 
dry veg 

fom rr 

green 

veg 

1064 

(APD) 

.02 

.88 

1 

.25 

.36 

.67 

.005 

.006 

.012 

liagHB 

.15 

.64 

.5 

.15 

.15 

1 

.004 

.004 

.003 

532 

(APD) 

.40 

.64 

.5 

.15 

.15 

.1 

.010 

.010 

.006 


Table 1: Return Rate FOM as a function of wavelength for various surface types and 
detectors and a clear atmosphere. Note that, except for green vegetation, the natural 
advantages of the 1064 nm wavelength with regard to surface reflectance and atmospheric 
transmission are largely overcome by the significantly higher quantum efficiency of 
detectors at 532 nm. 

The Return Rate FOM is the wavelength-dependent portion of Eq. (4), i.e. 

FOM m = k rj q pT*r] sh (24) 
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where X is the wavelength, q, is the detector quantum efficiency, p is the surface reflectance, T 0 is 
the one-way atmospheric transmission at nadir, and 

Photo-Diode (APD). 

Similarly, the Contrast FOM is given by the wavelength-dependent portion of (7), 

VshZ 


i.e. 


FOM c = 


JJsh _ = 


(25) 


Nl( 


1 + — 

8 P 




Nl(AA)A> 


. A , • the fiker p\VHM bandwidth, the factor of A 2 accounts for the reduced beam 

• *-* T and panmaer 4 u,cludes 

.1* : i w rpflftrtance and atmospheric transmission. 


A(nm) 

N x ° 

AA 

ang 

T]sh 

c 

soil 

dry veg 

C 

green 
ve S 

FOMc 

soil 

FOMc 
dry veg 

FOMc 

green 

veg 

1 1 

1064 

.06 

8.0 

l 

.94 

.95 

.98 

2.0 

Z.U 

>1 £ 

1. I 
A 0 

532 

/pirn 

.20 

1.5 

1 

.5 

.68 

.68 

.59 

4.6 

4.0 

H-.U 

a n 

\rivi i / 

532 

(APD) 

.20 

1.5 

.5 

.68 

.68 

.59 

4.6 

4.6 

4.U 


Table 2- Contrast FOM as a function of wavelength for various surface types an 
Note that the natural advantages of the 1064 nm wavelength with regard to surface 
re^ctance^and^olw background are more than overcome by the improved beam 
divergence and narrower spectral filters available at 532 nm. 

ildmete^orbiting 

planets «her than Earth, 532 ran would appear to be the wavelength of choice, 

2.8 EFFECTS OF DETECTOR AND/OR RECEIVER “DEAD-TIME” 

Following detection, teal detectora/timing fc£e«ed. 

within a few nsec), any dead time which is ong P )aser fire ( see Appendix B). 

be much reduced re,atlve to ** predicted by Eq 
(7a) for mean signal strengths substantially greater than unity. 
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Different photon-counting detectors have vastly different dead time characteristics. Avalanche 
Photo-Diodes (APD's) operated in the so-called Geiger-mode (also known as a Single 
Photoelectron Avalanche Diode or SPAD) trigger off a single photon and, following avalanche, 
must be restored via active or passive quenching techniques [Cova et al, 1996], The saturated 
output of a SPAD is suitably shaped for direct input into an event timer without an intervening 
timing discriminator, and hence the effective dead time is the longer of the timing receiver dead 
time or the detector quenching time which is typically 10 to 40 nsec depending on APD thickness 
[Cova et al, 1996], Since a photon entering a MicroChannel Plate Photo-Multiplier Tube 
(MCP/PMT) only depletes the charge in a few microchannels in the proximity of the cathode 
entry point, other photons (signal or noise) can be detected and amplified by the vast number of 
remaining unsaturated channels. Similarly, conventional dynode chain PMT’s usually require the 
input of many photons to saturate the output. Hence, most photomultipliers have virtually no 
internal dead-time, and the effective dead-time will generally be governed by the burst-mode 
characteristics of the timing receiver. 

It is demonstrated in Appendix B, that if the detector/receiver deadtime, r* is longer than the 
return waveform from the surface but shorter than the range bin duration (x b > x d , >c s ), the 
expression for the minimum signal cell contrast becomes 


c„ = i+(i 





(26) 


In (26) both the signal and noise counts per range bin, n s and n b , depend linearly on the receive 
aperture, A R . In the small aperture (signal) limit, the contrast approaches a maximum value 


lim C d = 1 + — 

^R 


(27) 


equal to the zero dead time result in (7a). However, in the large aperture limit, we obtain a 
contrast less than unity, i.e. 


lim C, = 1 — 


A* 


(28) 


suggesting that, if we choose too large a telescope, the signal will be indistinguishable from the 
noise regardless of the transmitted energy. 


In designing a practical instrument (see Section 3), we will often choose n 3 > 1 in order to get a 
high probability of detection per laser fire in single photon mode, making the exponential m ( ) 

small compared to unity. Thus, in the high signal limit, (26) reduces to 


lim = 1 - — + — 

n b 


(29) 
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As anticipated, comparing (29) with (27) demonstrates that the receiver effectively detects at 
most one photon in the signal swarm as a consequence of having a detector/receiver dead time 
longer than the return waveform. Any additional signal photons impinging on the detector are 
therefore wasted. Eq.(29) further indicates the importance of maintaining a background noise 
count per range bin, «*, significantly less than unity 

Solving for n h in (27) leads to an expression for the receiver area that maintains the desired 
contrast, i.e. 




P 


1 ) 


+ T, 


(30) 


where P is the maximum background noise rate per unit area per unit time defined in (5). 

Thus, reducing the receiver spectral bandwidth and spatial FOV to the maximum extent possible 
is an imperative first step in designing an efficient photon-counting altimeter. This avoids 
unnecessarily constraining the receive aperture and placing more of a burden on the laser 
transmitter to achieve the required power-aperture product for a given mission. As we will 
demonstrate in Section 4, the fundamental instrument parameters (transmitter energy, receive 
aperture, detector QE, etc.) must be chosen based on achieving the required sampling rate from a 
surface having a nominal worst case reflectance and slope. Over higher reflectance terrain or 
more modest slopes, the sampling rate will approach a maximum determined by the laser fire 
raters. We will demonstrate in later analyses that, over highly reflective terrain (e g., snow or 
ice), the increased solar flux from the surface can cause the noise count per range bin to rise to 
unacceptable levels, leading to a severe reduction in signal cell contrast (due to receiver dead time 
effects) in spite of the correspondingly higher signal levels. However, this can be accommodated 
by the inclusion of an actively controlled attenuator in the receiver to maintain system contrast 
and mimic the returns from a lower reflectivity surface. 


2.9 DETECTOR ARRAYS AND SEGMENTED ANODE PHOTOMULTIPLIERS 

Let us now assume that the terrain within the receiver FOV is imaged onto N p pixels of a photon- 
counting "array" detector and that each of the N p pixels is fed into a separate timing channel as in 
the airborne system to be described in Section 3 of this paper. Possible detector configurations 
include an array of Geiger-mode APD's, perhaps with a microlens array in front to reduce the 
amount of dead space in the detector plane, or alternatively by segmenting the anode of a 
conventional dynode-chain or microchannel plate photomultiplier [Degnan and McGarry,1997], 
One advantage is that each pixel is now monitoring the retums/solar noise from a much smaller 
area on the ground, i.e. A p = ApovANp . Furthermore, multiple signal counts can now be recorded 
since the detector and range receiver dead times in different timing channels are decoupled. Even 
if the detector elements are fed to a common event timer, a fixed time delay between channels can 
be introduced to overcome timer dead times and later removed from the range measurement in 
software while still preserving the identity of the active pixel [Degnan and McGarry, 1997], 

Since the receiver response is generally flat over the FOV, the noise counts are expected to be 
fairly uniformly distributed over the N p pixels or anode segments. The transmitted laser beam, on 
the other hand, typically illuminates a somewhat smaller region within the receiver FOV and is 
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most heavily concentrated near the center of the beam. The signal counts therefore will fall 
nominally on Ni pixels where Ni < N p . Thus, in the small aperture limit, the cell contrast in a 
given signal pixel channel (or the collective sum of all the active signal channels) is improved 
relative to the single element detector by roughly a factor of N/Ni as can be seen from (27). Of 
course, if aU of the independent channels within the receiver FOV are "ORed" together in the 
receiver prior to computing the cell counts, the cell contrast is unchanged, relative to the single 
element detector, in the small signal/small aperture limit. 

In the high signal limit (n, »Ni) , however, we see from (29) that, relative to the single element 
detector, the cell contrast for an individual signal channel is increased by roughly a factor of N p 
due to the fact that the solar noise counts are now shared equally by die N p detector elements 
within the receiver FOV but each of the Ni pixels viewing the ground laser spot is still capable of 
detecting only one signal photon due to receiver dead time effects. Furthermore, up to Ni 
individual photon TOF's can be recorded in a single laser fire interval, resulting in a single shot 
quasi-imaging capability provided, of course, that there are sufficient signal photons impinging 
on the detector. Thus, the instrument experiences an overall increase in contrast given by 


lim C d 


= 1 — 


• + 






(31) 


This suggests that an instrument providing approximately one photon per active pixel from the 
lowest reflectance surface of interest will, from a power efficiency standpoint, be near-optimum 
for rapid topographic mapping. We will demonstrate this for a Mars topographic mapper in 
Section 4. 

Solving for n b in (31) demonstrates that use of a pixellated detector in the high signal limit allows 
the receive aperture to be increased by roughly a factor of Ni relative to the single element case in 
(30) while maintaining a constant contrast, i.e. 


A r 


N, 


P 


(C d ~ 1) 


\-e 


+ t a 


(32) 


A further important advantage of the array detector is that the broadening of the return waveform 
(Poisson generating function) per pixel by beam divergence or surface slope effects is reduced by 
roughly a factor of 1/V N,. This, combined with the ability to more accurately define the source of 
the returning photon within the ground spot, results in a more accurate "quasi-point-to-point" 
range measurement. 

While they offer high angular resolving power for locating the photon source, imaging 
photomultipliers of the type discussed in the Introduction cannot be viewed as pixellated or 
segmented detectors. They have a single timing output, and the angular source of the photon in 
the far field is determined by additional circuitry. Thus, these devices must be treated as single 
element detectors, i.e. Ni = 1. 
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3. THE NASA AIRBORNE MULTIKILOHERTZ MICROLASER ALTIMETER 


3.1 PROGRAM BACKGROUND AND GOALS 

Under the auspices of NASA’s Instrument Incubator Program (HP), the Goddard Space Flight 
Center is developing an airborne version of the microaltimeter as a stepping ; stone to future space 
missions Instrument development was initiated in December 1998 and the first calibration flig 
are scheduled for January 2001. The goals of the IIP program are to: 

• Develop the necessary technologies and the operational and analytical software. 

• Confirm the validity of the mathematical modeling and surface data extraction 

algorithms. . n , 

• Collect high spatial resolution topographic data over a variety of surfaces (land, ice, 
water biomass, clouds) from a high altitude aircraft under both night and day conditions 
and evaluate its scientific usefulness relative to data from conventional altimeters. 

• Demonstrate the technical and economic advantages of the microaltimeter concept for 
future Earth and interplanetary missions. 

In order to improve the compatibility and portability of the instrument between different host 
aircraft, the microaltimeter has been designed to attach to a standard Lyca camera mo , 
commonly used in aerial surveying and photogrammetnc missions, and to operate at t^ica 
aircraft cruise altitudes between 8 and 12 km (25,000 to 40,000 ft). The key characteristics of the 
IIP prototype instrument are listed in Table 3. 


Transmitted Laser Pulse Energy 

10 mJ 

Operating Wavelength 

532 nm 

Laser Fire Rate 

10 kHz 

Transmitter Full Beam Divergence 

50 |xrad (1 

FWHM) 


Aircraft Altitude 

8 to 12 km (cruise) 

Ground Spot Diameter 

40 to 60 cm (depending on altitude) 

Receive Telescope Diameter 

1 4 cm (87% aperture utilization) 

Dftprtor OR 

12% (1,4, or 16 pixels) 

Detection Threshold 

1 pe (photon-counting) 

^nert ral Filter 

0.3 nm FWHM 

Receiver FOV 

| 140 prad 

(FWHM 


Panfry Resolution 

5 cm (detector limited) 


20 Hz 

bean rrequency 


Circular Scan (rotating optical wedge) 

300 to SOU m raai 

— 5 

US vuepenuiiiK uu aimuuv/ 

Effective Power-Aperture Product 

0.0013 W-m 


Table 3: Nominal Operatic Parameters of the NASA IIP Multikilohertz Mierolaser 
Altimeter 

3.2 REPETITION RATE CONSIDERATIONS 

The airborne microaltimeter is designed to operate at 10 kHz single photon sampling rates from 
aircraft cruise altitudes between 8 and 12 km (25,000 to 40,000 ft). While other a.rbomelaser 
altimeters [16] have operated successfully at these very high rates and have provided superb 
topographic nlaps, they typically use laser energies at least an order of magnitude higher and 
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operate at very low altitudes (1 to 2 km) where special FAA flight clearances are required Thus, 
our choice of a 10 kHz fire rate was driven by several factors: 

1 . We wanted to match the return rates of the best conventional altimeters but with a lower 
power-aperture product instrument operating at higher altitudes in order to demonstrate 

the technical and economic advantages of the microaltimeter approach. 

2 The 10 kHz rate allows system demonstrations to be carried out at aircraft cruise altitudes 
as high as 15 km (50,000 ft) without having to contend with multiple pulses in flight; this 
simplifies the receiver design. 

3. It appeared to be a comfortable rate for modem data acquisition systems and would be 
compatible with the projected near term capabilities of space-qualified CPU's or 
Programmable Array Logic (PAL) processors. 

The single pulse energy of 10 pj and 14 cm telescope aperture provide a mean signal strength of 
about 1 .5 pe per laser fire off a surface having a low reflectance on the order of 0. 1 and a mean 
slope of approximately 45°. This results in a 76% probability of detection for the aforementioned 
reference slope and approaches 96% sampling for level terrain with a nominal soil reflectance of 
.15 so that we can achieve the nominal 10 kHz surface sampling rate. As illustrated in Figure 1, a 
somewhat more efficient approach for the same power-aperture product would be to operate a 1 
pj laser at 100 kHz (for a return rate of about 0.1 pe per laser fire) and accept a roughly 90% data 
dropout, but this would have unnecessarily complicated the aircraft receiver and data acquisition 
by forcing it to track multiple pulses in flight at the much higher 100 kHz rate. The latter issue 
will be addressed in our later discussion of multi-kHz spacebome systems (Section 4) where the 
tracking of multiple pulses in flight is inevitable. 

3.3 INSTRUMENT OVERVIEW 

A block diagram of the NASA IIP instrument is presented in Figure 5. The system computer 
enables the firing of the passively Q-switched, frequency-doubled Nd:YAG microlaser 
transmitter and controls the repetition rate by adjusting the DC current to a 2.5 W CW laser diode 
pump array. The diode pump light travels to the laser via an optical fiber having an emitting 
diameter of 100 mm and a Numerical Aperture (NA) of 0.22 aid is refocused into the few mm 
long monolithic microchip by a short focal length lens. The microchip consists of a few mm long 
segment of NdrYAG (the gain medium) diffusion-bonded to, or grown onto, a thin crystal of 
Cr 4+ -doped YAG, which acts as a passive saturable absorber Q-switch and produces sub- 
nanosecond pulses. The matching indices of refraction and diffused nature of the interface results 
in low internal optical loss. The monolithic laser resonator is formed by polishing opposite sides 
of the microchip optically flat and parallel and applying reflective coatings [Zayhowski, 1998], 
The resulting microlaser, which is only a few mm in length, produces a train of low energy (~ 10 
pj), sub-nanosecond pulses at multi-kHz rates with no active switching components and, because 
of the monolithic structure, can never go out of alignment. 

A small pick-off mirror injects the transmitter beam, following some initial divergence, into the 
central 5 cm of the common transmit/receive telescope, which uses a 14 cm off-axis parabola as 
the telescope primary mirror. Optical leakage from the outgoing laser pulse is sampled by a fast 
photodiode, which time-tags the start pulse with a Berkeley Nucleonics Corporation (BNC) 
Time-to-Digital Converter (TDC) and starts the range gate generator (BNC Model B940). The 
TEMoo spatial mode quality (M 2 <1.3) of the end-pumped microlaser allows an extremely narrow 
transmit beam (~10 arcsec) to be generated with a correspondingly small receiver FOV, which in 
turn reduces the solar noise count rate. The laser spot diameter on the ground varies between 40 
and 60 cm, depending linearly on aircraft altitude. The receiver utilizes the remaining 87% of the 
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available aperture that is not blocked by the injection mirror to gather photons reflected fromthe 
surface. Stray light rejection within the receive telescope and interface optics severely restricts the 
number of off-axis photons reaching the stop detector. Internal gating of the detector and range 
receiver also helps prevent saturation of the detector resulting from laser backscattenng within 
the instrument or aircraft. 


The returning photons pass through spectral (0.3 nm) and spatial (-140 prad) filters and are 
imaged onto a 2x2 element Hamamatsu Model R5900U metal channel dynode Photo-Multiplier 
Tube (PMT) Because the primary timing unit, the BNC Model B945 TDC, can accommodate up 
to 16 stop events (one in each of 16 channels), the receiver design permits 4 stop events to occur 
within the range gate for each of the 4 detector elements. Provisions have also been made in toe 
design to incorporate either a single element Model R7400U (1 channel, 16 stops) or a 4x4 Model 
H6568 PMT (16 channels, 1 stop per pixel) as alternate detectors. Use of pixellated 
photomultipliers permit a crude, but inexpensive, demonstration of the potential advantages of 
"quasi-imaging" or "point-to-point ranging". For example, at 12 km altitude where the ground 
beam diameter is about 60 cm, the 4x4 array allows the source of the received photon to be 
resolved within a roughly 15cm x 15 cm cell on the ground. In a 300 km Mars orbit, the r ° u 8 h| y 
15 m diameter beam would be dissected into 4 x 4 meter cells by a 4x4 array leading to a higher 
resolution ground "image". Of course, improved spatial resolution could be achieved either 
through higher pixellation or the use of imaging photomultipliers as discussed previously in toe 
Introduction. It is worthwhile to note that, to the extent that the signal and noise generation are 
both reasonably uniform within the receiver FOV, splitting the receiver FOV into individual 
pixels has no impact on the instrument contrast since the expected signal count and noise count 
per pixel are reduced equally. 

Each detector pixel output is input to a high speed amplifier/constant fraction discriminator 
(AMP/CFD) module (Phillips Scientific Model 6908). The CFD provides a fast tuning pulse to 
both a "coarse" and a "fine" range receiver. The function of toe coarse receiver is to 
simultaneously capture all of the "soft" scattering surface returns (e.g. clouds, boundary or fog 
layers) as well as the ground return and to provide an initial estimate of where the ground and 
near-ground scattered (e g. tree canopies and sub-canopies or man-made structures) are located 
so that the range gate on the fine receiver can be substantially reduced and optimally placed m 
time to concentrate on these surface features. 

The "coarse" receiver is an EG&G Ortec Turbo- Multi-Channel Scalar (MCS) which bins all of 
the photons received over virtually the entire fire interval and creates a lustogram ° f P h ?^° n ^ p 
events. The best coarse range resolution of 75 cm is determined by the minimum MCS bin size of 
5 nsec Since this level of range accuracy is more than adequate for clouds or other 
meteorological layers and there is similarly no need for high spatial resolution m toe transverse 
dimension, the "stop" outputs of all toe pixels are ORed together at toe input to toe MCS. 

The "fine" range receiver time-of-flight (TOF) measurement is determined by differencing toe 
time tags of toe laser fire ("start") event and the photon return ("stop") events in the TDC. 
Although toe TDC has a timing resolution of 50 psec (7.5 mm), toe range resolution of toe 
receiver is limited by toe detector impulse response to about 5 cm RMS Furthermore, since toe 
TDC is capable of measuring only one stop event per channel, toe serial stop pulses (up to a 
maximum of four) from each of toe four pixels must first be separated into p arallel taesby a 
channel Serial-to-Parallel Converter (SPC). The SPC output fens out to the 16 channels of the 
TDC, thereby preserving pixel identification. 
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In order of decreasing duration, individual gating pulses are applied to the PMT, the MCS, and 
the TDC respectively by the multi-channel range gate generator as in Figure 4. In initial flights of 
the prototype instrument, all photon “times-of-arrival” within the range gate, relative to the start 
pulse, will be measured by the "fine" range receiver and recorded on the system hard disk for 
later analysis. However, our data extraction algorithms will be running in parallel to validate them 
in preparation for later flights, where software algorithms will identify the signal photons in 
flight and strip away most of the noise photons in near real-time. 

There are three computers used in the flight system: the Ranging System Operator Interface 
(RGUI), the Ranging System Data Collector (RDATA), and the Navigational and Camera 
Control (NAV). A common timing reference is input to all three computers in the form of the 
current laser shot counter, which allows post-flight coordination of data recorded on separate 
computers. 

The Ranging System Operator Interface is a 500 MHz Pentium III processor running Windows 
NT 4.0. The system has 128 Mbytes of memory and two 13 Gbyte IDE Ultra ATA hard drives. 
One of the drives is used for system and program files; the other is mounted in a receiving frame 
for easy removal and will be used to record the flight data. The operator display is an LCD 20 
rack mounted monitor with a pixel resolution of 1280 x 1024. The RGUI computer is responsible 
for operator interface, data logging, signal processing, and control of the EG&G Turbo-MCS 
Histogrammer. In addition this computer picks up a small data set from the NAV computer at 1 
Hz using internet file sharing, and it reads 720 Kbytes of shared memory data from the RDATA 
computer at 1 Hz. Signal processing in RGUI provides visual feedback to the operator and can be 
used to provide control of the range window. All of the data is logged to the removable hard disk 
once per second. The amount of data recorded is approximately 3 Gigabytes per hour. 

The basis of the RGUI software is NTGSE, which is a software package originally developed at 
GSFC for other missions. NTGSE has been modified for our application and has allowed a 
quicker software development cycle. 

The Ranging Data Collector (RDATA) is a 233 Mhz Pentium processor running DOS 6.22. The 
processor is an embedded module in a National Instruments VXI crate, which also contains a 
MXI interface to the RGUI computer (with 32Mbytes of shared memory), a digital I/O module 
(National Instruments DIO- 128), and two Berkeley Nucleonics modules for ranging interface: the 
B940 digital delay generator (DDG) and the B945 1 6-channel time to digital converter (TDC). 

The embedded CPU contains 32 Mbytes of memory, a 4 Gbyte hard disk, as well as ethemet, 
SCSI, serial and parallel interfaces. RDATA uses an interrupt from the TDC module at the end of 
each range gate to trigger the reading of the laser fire time, die 16 range returns, the scan angle, 
and the shot counter. Interrupts are generated at the laser fire rate of 10 kHz. Data is placed in a 
shared memory buffer for the RGUI computer to pick up when it can. The shared memory buffer 
is a circular buffer and currently holds 5 seconds worth of data. Commands from the RGUI to the 
RDATA computer are also sent across shared memory. These include the range delay and the size 
of the range window, and are picked up by RDATA at 20 Hz. 

The instrument will be flown with and without a simple scanner in order to demonstrate both 
contiguous 2D linear mapping of the underlying terrain as well as the 3D mapping of large 
swaths, extending over several hundred meters in the transverse direction, during a single pass of 
the aircraft. Because of the small size of the microaltimeter transmit/receive optics, an 
inexpensive optical wedge roughly 15 cm in diameter can be rotated at 20 Hz to superimpose a 
circular scan on the linear flight path motion. The rotating wedge, which deviates the beam from 
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nadir by approximately 1 degree, simultaneously deflects both the transmitter beam and narrow 
receiver FOV while maintaining a common boresight 

3.4 POST-FLIGHT DATA PROCESSING 

In order to generate an accurate topographic map, it is essential for any laser altimeter that we 
know both the instantaneous position and attitude of the instrument on each laser fire in addition 
to the pulse time-of-flight and any fixed timing biases within the instrument. The latter bias can 
be determined by ranging to a target whose distance from the instrument reference point is well 
known. For the quasi-imaging and scanning microaltimeter, one must also make corrections 
based on the instantaneous off-nadir angle of the scanner and the X-Y coordinate (or pixel) which 
recorded the photon event 

Multiple GPS geodetic receivers will be used (one in the aircraft and one or more on the ground) 
to provide post-flight dynamic differential positioning at the few decimeter level. The fixed offset 
vector between the GPS antenna phase center and the altimeter reference point in the aircraft 
reference frame will be measured and corrected by three-axis attitude data to accurately locate the 
altimeter reference point in the Earth reference frame during flight. Instantaneous instrument 
attitude will be provided by a three-axis fiber-optic gyro mounted directly on the altimeter optical 
bench and calibrated pre-flight. To counter or monitor gyro drift errors, independent updates of 
attitude will be obtained in-flight by (1) a digital compass (for heading) supplemented by two 
collocated and orthogonal bench-mounted inclinometers and (2) repeatedly over-flying a set of 
four "ground stars", whose relative positions have been precisely located at the few mm level 
using geodetic GPS receivers. These "ground stars" will be located near the center of the 
interrogated region to accommodate frequent overflights and will be arranged in a square or 
diamond pattern a couple of km on a side. Thus, when the aircraft is near the center of the pattern, 
they can all be viewed simultaneously by the wide FOV camera attached to the altimeter 
instrument, providing adequate angular resolution for updating the aircraft attitude at the few 
arcsecond level. 

One candidate approach to generating the ground stars is to use small battery-powered Light 
Emitting Diode (LED) arrays, which are pulsed by the lpps output of a GPS Timing Receiver . 
The onboard digital camera, boresighted with the range receiver and synchronized with the wide- 
angle ground star emission by the lpps output of the onboard GPS receiver, then records the 
spectrally narrow images of the ground stars through a bandpass filter. The deviation of the 
ground stars from their predicted positions in the camera FOV, based on the few decimeter 
accuracy aircraft positioning, allows the instrument attitude to be updated post-flight with few 
arcsecond accuracy. 


The instantaneous angular position of the scanning wedge is obtained by interpolating between 
multiple sensors which monitor the passage of a particular point on the wedge by producing an 
output pulse which is time-tagged and recorded in the data file along with the ranging data. This 
was viewed as an adequate and less expensive and data intensive approach than using encoders 
due to the fairly constant wedge velocity between monitoring points. As mentioned previously, 
the detector pixel producing a given photon event is recorded by the "fine" range receiver and can 
be used to correct the measured elevation for the transverse location of the photon source within 
the ground spot. 

To summarize, a "start" pulse and potentially one or more "stop" pulses (noise and/or signal 
counts) will be recorded, for each laser fire, by the coarse and fine receivers, both of which can be 
viewed as correlation range receivers but with vastly different range gates and resolutions. The 
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post-detection Poisson filtering techniques described in feft^ 

roll and vaw in order to determine the terrain heights. Additional corrections i * 

measured terrain heights in the transverse dimension. 

SSSESsSraSSfar 

differences and an overall measure of terrain difference 

Follow-on science flights will investigate instrument performance ^s^ta ^ 

ss‘x«.=SS=- 

2^^^£Sss=5s=sss* 

interpretation of the high repetition rate data, own frequent and large elevation 

instrument and date processing algontems to of beach areas will 

SESES* fcrto™ia“n and wave heigh* and for performing shallow 

water bathymetry at single photon levels. 

M SIMULATION OP THE AIRBORNE MICROALTIMETER PERFORMANCE AT 12 KM ALTHUDE 

Our analysis of the airborne microahimeter is summarized by the six plots in 5“* 

?£££ three curves, labeled by tee assumed ? 

slope of 45 degrees was assumed J^S m Appendix A 
in a single shot return waveform of about 4 ns mos! terrain within a 4 

indicate that 2 to 4 noise counts per laser fire are e>?e«dto °™'^ r Over highly 

microsecond gate which is well within the nsec 

reflective clouds or snow, this can grow to 11 < A nominal 4 
was estimated for the aircraft hardware so a "fine" receiver range window of 

microsecond timing gate was assumed, w ic co sp , - mapping of the tallest forests 

600 m. the latter value would acconmo^e detailed to wtei i>W « (the 

detector. 

Figure 6a plots the maximum surface aid the 

altimeter operating at the marnurn ^ ls ^ a s ^ d from ^ figure ^ minimum power-aperture 

pr^^quiredto achieve 10 kHz surface l 

^ sf detection on each laser fire with a nominal 


22 



mean signal strength of l pe per laser fire (see Figure 1) and to avoid the need to operate the 
system at 100 kHz rates as described previously in section 3 .2. 

Figure 6b plots the minimum cell contrast as a function of the number of laser fires per frame 
(frame interval). The cell contrast decreases with increasing frame interval because ofthe need to 
widen the range bin according to Eq. (23a). Over the range of 1 to 100 laser fires per frame, 
which is common to Figures 5(b) though 5(f), the frame interval increases from Orf to 10 msec, 
the range bin width varies from 4 to 10 nsec, the corresponding number of range bins decreases 
from 1000 to about 400, and the ground distance covered by the aircraft at a velocity of 10 sec 

varies from 1 cm to 1 m. 


Figure 6c plots the optimum frame threshold as a function ofthe number of laser fires per frame 
as computed from (19) and rounded off to the next highest integer. 

In Figures 6d through f, we plot the Differential Cell Count (normalized to M), the probability of 
correctly identifying the signal cell, and the mean number of falsely identified cells per frame 
respectively as a function ofthe frame interval. We note that the Normalized Differential Cell 
Count approaches the ideal value of unity for 15 laser fires per frame or a frame interval o 
msec. As expected, the probability of detecting the signal cell also approaches unity and the mean 
number of falsely acquired noise cells per frame fells well below unity for these longer frame 
durations. However, one should choose the minimum frame duration that gives acceptable results 
since ultimately the range bin can grow to encompass a large fraction ofthe range gate w ic 
limits the effectiveness ofthe selection algorithm as a real-time noise filter. 

It should be noted that these analyses assume a worst case solar noise background as calculated m 
Section 2.4. The frame threshold and frame interval can be initially set based on a pnon 
expectations and then recomputed in flight using (19) and actual measurements ofthe mean 
signal and noise counts per range cell. A smaller mean noise count will increase the contra J 
which in turn reduces the threshold and the required frame interval. Adaptively 'narrowing the 
range gate over barren land or low vegetation to reduce the value of N to in (19) has a similar 

effect. 

Figure 7 shows the results of a Monte-Carlo simulation ofthe microaltimeter operating in non- 
scanning mode over a rapidly changing terrain devoid of any ground vegetation. The simulated 
surface which includes a 50 meter cliff, is easily seen against the solar background noise in the 
coarse map of Figure 7a. In spite ofthe steep cliff, the selection algorithm was able to rapidly 
identify and reacquire the surface returns from the top ofthe cliff. 

The expanded scale in Figure 7b shows the actual roughness for a segment of fee simulated 
surface in Figure 7a as well as fee data points selected by fee algorithm as surface returns (points 
surrounded by circles). The few noise counts falsely identified as signal in Figure ^ happened to 
fall into fee chosen signal cells and would be eliminated in iterative applications of fee selection 
algorithm if one substitutes fee bin size ofthe signal cell for fee range gate and subdivides it into 
still smaller range bins on each iteration. 
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4.0 SPACEBORNE ALTIMETERS FOR PLANETARY TOPOGRAPHIC MAPPING 

4.1 ADDITIONAL TECHNICAL CONSIDERATIONS FOR SPACEBORNE MISSIONS 

To simplify our discussion of spacebome microaltimeters, we will consider a planet other than 
EarT^ose surface topography is umquely complicated by the presence ofvegetaDon ^ man- 
made structures. We will also consider a non-scanning system with a single element drtector. In 
scaling the microaltimeter from airborne to spacebome platforms, we encounter the following 
technical issues. 

1) The required power-aperture product for the altimeter scales as R* for a given surface 
sampling rate, but, if the receiver dead time is long compared to the single shot surface 
return waveform, the receive aperture will necessarily be limited by the need to keep 

noise counts within acceptable limits. 

2) One must contend with multiple pulses in flight in order to achieve kHz or higher 

sampling rates from orbit. , . . „„„„ . - 

3) Higher spacecraft ground velocities imply a faster rate of growth in the required range bm 
with increasing frame interval and modest surface slope. In addition, the roundtnp transit 
time to the surface is on the order of milliseconds which contributes to the uncertainty m 
the predictive range window and further increases range bin duration. 

4) Prime power, size, and weight play a more significant role in spacebome instrumentation, 

especially for interplanetary missions. 

We will address each of these topics as we proceed with an analysis of a Mars topographic 
mission. 

4.2 ONE METHODOLOGY FOR DESIGNING A MULTI-KHZ MARS MICROALTIMETER 

Using our theoretical results from Section 2, we will now illustrate one methodology for 
designing an optimized multi-kHz microaltimeter to rapidly map the surface of Mars. The des gn 
process begins by: (1) selecting the minimum sampling rate one wishes to achieve over a nomrn 
worsTcasesurface (i.e, maximum median slope and mmimum refkmPmce) tPI— ^ 
instrument optical losses; (3) minimizing the norse background vta 8 

and utilizing a correspondingly low transmitter beam divergence; and (4) selecting a high 

quantum efficiency detector 

A typical altitude for a Martian satellite (e g., MOLA) is 300 kM with a «>^ s P OTdin8 i 5°“"^ rv 
velocity slightly greater than 3 km/sec. Thus, if we wanted to sample the worst case smfece every 
meter along-track, we would set/, • 3 kHz. We note from (3) that, u, commode the 

sampling rate,/,, is usually less than the laser fire rate,/,, since the probability of detecting a 
signal photon is less than unity unless «, »1 

Assuming a 50 prad foil beam divergence, which is easily achieved with a transmit aperture of 
one to two inches and the TEM00 mode quality of the microlaser source [Zayhowski, 1998], the 
beam diameter on the ground is 1 5 m. Thus, the mean spacmgbetweensubsequeirt^crts entire 
ground is about 7% of the beam diameter and the returns are therefore highly contiguous We will 
forther assume that the receiver FOV is twice as large (100 prad) as the transmit beam m order to 
accommodate some alignment and/or stability errors m the boresighting of the transmitter and 

receiver. 
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The median slope measured by MOLA, on 300 m beam centers along-track and averaged over 35 
km baselines, as analyzed by MIT rarely exceeded 5 degrees and the legend on their global map 
shows a maximum plotted value of about 23 degrees [Aharonson, 2000], which we will adopt as 
We will also assume that the nominal reflectance of soil on Mars is comparable to that on 
Earth, i.e. about 15% at 532 nm. Locally, the slopes and reflectances may vary about these values 
leading to both increases or decreases in the photon return rate from the surface, but our auto- 
tracking algorithm is capable of adjusting to these variances. 

Over high reflectance surfaces (if any), the mean signal photon counts will increase, and the 
microaltimeter sampling rate will approach a maximum determined by the laser fire rate. 

However, because of the receiver dead time effects discussed in Section 2.8, we must include an 
optical attenuator in our receiver design in order to reduce the mean noise count per cell and 
maintain the system contrast at an acceptable level by mimicking a lower reflectance surface. 
Without the attenuator, a sufficiently bright surface would be indistinguishable from the solar 
background. 

The GLAS project plans to use a space-qualified photon-counting detector with a quantum 
efficiency greater than 50% at 532 nm [Abshire et al, 2000]. The SPAD detector is actively 
quenched with a dead time of 50 nsec. The principal optical loss element in the receiver is the 
narrowband spatial filter. Spectral bandpass filters with good optical throughput (>50%) at 532 
nm are generally limited to a spectral bandwidth of AX — 3 A° FWHM. We will adopt this value 
for the purposes of this analysis and assume a total receiver optical throughput efficiencies of rj r - 
0.4. If necessary, spectral filtering can be further enhanced, via a thermally tuned etalon filter, to 
as little as 0.3 A 0 with an additional 15% loss in throughput as in the GLAS instrument [Abshire 
et al, 2000], 


Instrument C 

haracteristics 

Wavelength, X 

532 nm 

Detector Quantum Efficiency, 

0.5 

Detector Pixels, N„ 

1 

Receiver Optical Efficiency, rjr 

0.4 

Transmitter Full- Angle Divergence , 26, 

50 prad 

Receiver Full Angle FOV, 26 r 

100 prad 

Spectral Filter Bandwidth, AX 

3 A 0 

Receiver Dead-Time, Td 

50 nsec 

Noise count rate per unit receive area, (5 

5 x 10'/m 2 -sec 

Minimum sample rate, f s 

3kHz 

Spacecraft Characteristics 

Altitude, R 

300 km 

Ground Velocity, vf° 

3.13 km/sec 

Surface Characteristics 

Reflectance, p 

0.15 

Maximum Median Slope, cr max 

23° 

Atmospheric < 

Characteristics 

One-Way Transmission at Nadir 

0.90 


Table 4: Summary of parameters used in Mars micro altimeter analysis. 
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Substituting the aforementioned values, which are summarized in Table 4, into (4) and solving for 
the minimum power-aperture product yields a value P A mw = 0.017 W-m . For night only 
operations, one can arbitrarily trade off receive aperture and laser power to achieve the necessary 
product. However, we will now demonstrate that, in the presence of solar noise, this minimum 
power aperture product, which from (3) implies a mean signal count per laser fire on the order of 
0.1 pe, is not necessarily the optimum choice. In order to reduce the mean noise count per range 
bin well below 0. 1 pe and maintain adequate cell contrast, the receive aperture must be kept 
relatively small, placing a greater burden on the laser transmitter power to achieve the required 
power-aperture product. In all potential instrument configurations, the receive aperture 
requirements are fairly modest (10 to 20 cm), and one may wish to choose an instrument 
configuration which places the fewest demands on the laser. 

Since Mars lies a mean distance of 1.524 AU from the Sun, the spectral irradiance, impinging 
on the planet's atmosphere is reduced by a factor of 2.32 relative to Earth to a value 0.086 W/m - 
A 0 at 532 nm. Furthermore, since Mars' thin atmosphere is expected to be more highly 
transmissive than Earth's, the atmospheric scattering component of the noise is reduced relative to 
the surface component. For the purposes of our analysis, we will assume an atmospheric 
transmission of 90% . From (A8) and (All) in Appendix A, we calculate a worst case solar count 
rate per unit receive area of 4.2 x 10 7 /m 2 -sec from the Martian surface and 8.2 x 10 /m -sec^from 
the intervening atmosphere. Thus, the total worst case solar noise count rate is fi— 5.0 x 10 /m - 
sec with 84% due to surface scattering. If we assume a typical receiver dead-time, r d , of 50 nsec 
in (30), we obtain an upper limit on the receive aperture ofA max = l/{$ r d = 0.4 m 2 , corresponding 
to a maximum receiver diameter of 7 1 cm in the low contrast limit. 

We demonstrated in Section 2.5 that different combinations of contrast and frame time can 
provide roughly the same ability to distinguish signal from noise. The required signal count per 
frame, N s , increases with decreasing signal cell contrast as demonstrated in Figure 4. The frame 
duration is given by the formula r f =N/f s . The corresponding range bin can be computed from 
(23) with y = 1 .77 x 10' 5 . The results are summarized in the following table: 


Contrast, 

C 

Signal Counts per Frame, N s 
Ground Track per Frame (m) 

Frame Duration, 
T r (msec) 

Range Bin, 
x b (nsec) 

2 

80 

27 

515 ~1 

3 

50 

” 17 

338 

5 

35 

12 

249 

10 

25 

8 

190 


Table 5: Combinations of contrast and mean signal counts per frame which provide roughly 
the same ability to discriminate signal from noise. Because for f t — $ kHz we make on 
average one range measurement for every meter of ground travel, the ground track covered 
per frame is roughly equal to N, in meters. The corresponding frame durations and range 
bins for the postulated Mars mission are computed using various formulas cited in the text. 

We now let the mean signal count per laser fire, n s , vary between 0. 1 and 10 pe and compute in 
succession the following quantities: 

Single-shot probability of detection: P d = 1 - exp(-n i ) 

f 

Laser Fire Rate: f os = — L from (3) 

* A 


26 



Receiver Area: A r - 


1 


from (30) 


C-l 


-T b + T d 


Transmitted Single Pulse Energy: E, = 


CA r 


from (1) 


Transmitted Laser Power: P, = f gs E, 


Power-Aperture Product: P t A r 


n, fj t 

\-e n - C a 


The results are plotted in Figure 8 as a function of the mean signal count per laser fire over a 
range 0. 1 pe <n, < 10 pe. In, Figure 8a, the single shot probability of detection is greater than 
98% for n, > 4. Figure 8b illustrates the fact that, at low mean signal strengths on the order of 0. 1 
pe, the altimeter must be operated at 30 kHz in order to achieve the 3 kHz sampling rate and 
compensate for the low probability of detection. The required laser fire rate falls asymptotica y 
to 3 kHz as the probability of detection asymptotically approaches unity with an increased mean 
signal count. In (c) through (e), the four curves are labeled by the assumed value for the ccmtrast, 
C, with the corresponding value of the mean signal count per frame, N s , given by Table 5. In 
Figure 8c, we use (30) to demonstrate that a stronger mean signal count allows the receiver 
diameter to grow until the mean noise count per bin per fire gets close to the limiting value, n b - 
1/(C-1). This results from the fact that the system dead-time allows the detection of only one 
signal photoelectron per laser fire as discussed in Appendix B. Since the receiver area cannot 
grow further once this limit is reached, the transmitted pulse energy must increase to achieve still 
higher mean signal counts as in Figure 8d. Figure 8e provides the interesting and useful result that 
the required average laser power reaches a minimum for all contrasts at a mean signal count o 
about 1 .24 pe. As the desired contrast increases, however, the required laser power also increases, 
again due to dead-time effects which constrain the receive aperture. Finally, the Power-Aperture 
product is plotted in Figure 8f, and, as expected, all four contrast curves overlap perfectly. As 
expected from the theoretical results in Section 2. 1, the power-aperture product required to 
achieve a given sampling rate increases monotonically with increasing n s . 

Returning to the point of minimum transmitter power in Figure 8e, we note that a transmitted 
laser power of only one watt, coupled with a very modest 18 cm diameter receive telescope 
would allow us to do nominal 3 to 4.3 kHz ranging to the Martian surface from a 300 km orbit, 
provided we accept a minimum signal contrast of 2. The laser would operate at roughly 4.3 kHz 
with a single pulse transmitted energy of only 230 pJ. Choosing a higher signal contrast of 10 
would allow use of a smaller 10 cm telescope, with perhaps a correspondingly smaller and lighter 
scanning mechanism (if relevant), but it also would require a larger and heavier 3 watt laser (700 
uJ (a) 4.3 kHz) with approximately three times the onboard power requirements. As demonstrated 
in Figure 8e, the power-aperture products would be identical for the two systems (0 024 W-m ), 
but the lower contrast system would place less of the signal strength burden on the laser. This 
would normally make it the system of choice since, by conventional spacebome altimeter 
standards, the telescope size and weight is very modest in either case. The principal advantage of 
high contrast systems is that fewer noise events are retained by the autotracking algorithm m real 
time, but, as discussed in Section 2.5, in-flight post-processing of the data can easily remove 
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these. As a point of reference, the power-aperture product for the MOLA altimeter is almost 4 
times larger, i.e. PA M ola = 0.088 W-m : , yet the sampling rate is 300 times smaller (10 Hz vs 3 
kHz) for the nominal "worst case" surface 

To summarize, we argued in Section 2. 1 that operating the altimeter at low mean signal levels (n s 
< 0. 1) allowed us to maximize the sampling rate for a given power-aperture product. In the 
absence of background noise, we are free to arbitrarily choose the transmitted power and receive 
aperture which yield this product. In the presence of noise, however, the receive aperture must be 
severely limited in order to maintain the necessary signal contrast (i.e. n b < n s < 0. 1), which in 
turn places more of the signal burden on the laser transmitter. By operating at a higher signal 
level slightly greater than 1 pe (and hence a somewhat higher power-aperture product), we can 
increase the size of the receive aperture and minimize the required transmitter power. 

4.3 TIMING OF MULTIPLE PULSES IN FLIGHT 

In the 10 kHz airborne system discussed in Section 3, only one pulse is in flight at a given time 
for aircraft altitudes below 1 5 km. Thus, a simple event timer, triggered by the outgoing pulse, is 
sufficient for the range measurement. For a 4 kHz microaltimeter operating from a circular orbit 
300 km above Mars, however, there are approximately eight pulses simultaneously in flight, 
separated by roughly 75 km. 

The problem of designing a correlation range receiver (CRR) which handles multiple pulses in 
flight has been addressed previously in connection with NASA's developmental photon counting 
SLR2000 satellite laser ranging system, designed to operate with segmented anode 
photomultipliers at a 2 kHz laser fire rate. An early version of the receiver [Degnan and McGarry, 
1997] had a demonstrated RMS ranging resolution of 4 mm and was later improved to about 1 
nun [Degnan, 1998], Since such high range resolution is unwarranted for an altimeter designed to 
have few cm resolution, one can consider simpler and more power efficient alternatives for the 
timing of photon "events". 

The range receiver is designed around a "sub-clock" which produces timing pulses at the chosen 
laser fire rate, f QS , and is derived through frequency division from the high rate output (typically 
10 MHz) of the spacecraft clock/frequency standard. This "sub-clock" pulse can then be used to 
trigger both the microlaser transmitter pump diodes and a timing receiver, which is operated 
modulo the laser fire interval (~ 250 psec for our 4 kHz Mars example). The pump diodes 
produce a rectangular optical pulse which pumps the Q-switched microlaser until laser threshold 
is reached, at which point an intense subnanosecond pulse is rapidly produced. The timing jitter 
between the laser pump trigger and the emission of the laser pulse can be large compared to the 
desired few hundred picosecond timing resolution. Thus, to achieve the necessary range accuracy, 
the temporal positions of the outgoing pulses relative to the sub-clock "tick" must be measured on 
each fire. The return photon from the Martian surface arrives 8 fire cycles later, and its time of 
arrival is again measured relative to the sub-clock timing pulse. By properly pairing the "start" 
and "stop" events, accurately measuring their temporal positions relative to the sub-clock "ticks", 
and adding the appropriate number of intervening sub-clock cycles, an accurate time-of-flight 
measurement can be obtained [Degnan and McGarry, 1997]. Thus, each photon event can be 
labeled by: (1) the epoch time second on the spacecraft clock preceding the measurement; (2) the 
laser fire interval within the following second (0 - 4000 for a 250 psec interval) during which the 
photons departed the instrument ("start event") or were received ("stop event"); and (3) the time 
interval between the sub-clock pulse and the photon event as measured by the modulo timer. 
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Because of the large 75 km spacing between pulses, the number of intervening clock cycles can 
be discerned from a priori knowledge of the satellite altitude obtained from navigation data. 
Furthermore, at 4 kHz, the 250 psec fire interval allows a range depth up to 37.5 km to be probed, 
which is long enough to easily accommodate the entire topographic range of the inner planets. 

For example, we note that the highest point on Earth, Mt. Everest in Nepal, lies roughly 8.8 km 
above mean sea level. Thus, all surface and canopy returns (plus most low lying clouds of interest 
for radiation transfer studies) could be accommodated within a fixed range gate on the order of 60 
psec if one were to assume a fixed receiver range gate, a highly circular orbit, no terrain below 
sea level, and a truly spherical planet. In reality, because of the planet’s oblateness, an altimeter in 
circular orbit about the Earth would require a slightly larger temporal window on the order of 75 
psec to accommodate the global surface returns. Still, this global variation in surface topography 
represents only 33% of the 250 microsecond period between pulses. Extending the range window 
to span either the entire 20 km height of the Earth’s atmosphere to include high cloud returns for 
radiation balance measurements or the highest feature on Mars, Olympus Mons at 2 1 km [Frey, 
2000], causes a global spread up to 146 psec or 55% of the laser fire period. Furthermore, since 
the laser fire is triggered by the sub-clock ticks and the jitter in the pulse output is also a relatively 
small fraction of the fire period, the start and step pulses can be temporally isolated from each 
other within the laser fire period via a fixed range gate generator and input to a single modulo 
timer to conserve prime power. 

At laser fire rates greater than about 7.5 kHz, the range return intervals from different laser fires 
will begin to overlap in time at the receiver for probed range depths on the order of 20 km, 
making data interpretation more difficult. To avoid this complication, higher repetition rate 
systems would need to reduce the range depth and concentrate on surface and/or near-surface 
returns (e.g. tree canopies and low lying cumulus clouds. 

Finally, if the 250 psec laser fire interval for our Mars Microaltimeter is broken up into 500 range 
bins, as assumed in Figure 4, the resulting range bin duration of 500 nsec is a good match for the 
C = 2 case in Table 5, which required the least transmitter power and hence was the instrument 
configuration of choice. 

4.4 A POWER EFFICIENT ALL-DIGITAL TIMER 

As mentioned in the Introduction, conventional laser altimeters often make use of high speed time 
digitizers which record the waveform of the multi-photon returns. When compared to digital 
circuits, these analog devices tend to consume large amounts of prime power as and require the 
use of complex (and often fallible) algorithms to compute the range to complex surfaces. On the 
other hand, computing the time-of-flight of a single photon is a relatively unambiguous process 
and, combined with the quasi-imaging capabilities of pixellated detectors, permits greater spatial 
resolution in both the horizontal and vertical directions. Furthermore, the modest ranging 
resolution (few ems) required by the altimeter allows one to consider a compact and power- 
efficient all-digital design for the modulo timer [Degnan et al, 1999], 

A block diagram of the All-Digital Event Timer, with an integrated 6x6 micro-APD array, is 
shown in Figure 9. A 4 GHz clock is derived from a voltage controlled crystal oscillator and is 
phase-locked to the 10 MHz frequency standard. A 21 bit synchronous counter is used to time tag 
photon events from the 36 element array. An input signal conditioner is used to "latch" the 
incoming photon events to the nearest 4 GHz clock pulse, providing a timing resolution of 250 
psec. The latch stores the counter data to time the event. An array location data generator 
identifies which pixel of the array has been excited by the incoming photon. The pixel location 
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information and the high precision timing information are then fed to a 64 bit FIFO (First In Fi 
out) buffer The resulting FIFO output is a 64-bit word that contains the time of the occurrence ot 
the event and the location of the excited pixel. We estimate that low voltage (<3 V) integrated 
circuits (ICs) and differential signal (LVDS) drivers and receivers can reduce power 
consumption of the timer to less than 6 watts in a very small package occupying about 5 cubic 

inches. 


5.0 SUMMARY 

We have described a photon-counting laser altimeter whose advantages include: (1) a more 
efficient use of available laser photons providing up to two orders of magnitude greater surface 
sampling rates for a given power-aperture product; (2) a simultaneous two order o magn e 
reduction in the volume, cost and weight of the telescope system; (3) the unique ability to 
spatially resolve the source of the surface return in a photon counting mode through die use ot 
pixellated or imaging detectors; and (4) improved vertical and transverse spatial resolution 
resulting from both (1) and (3). Furthermore, because of the significantly lower laser pulse 
energies, the microaltimeter is inherently more eyesafe to observers on the ground and less prone 
to the internal optical damage, which can terminate a space mission prematurely. 
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APPENDIX A: ALTIMETER NOISE MODEL 


A*1 SOLAR SCATTERING FROM THE ATMOSPHERE 

Because the effective height of the atmosphere above sea level (h, ~20 Km) is small ^lativeto 
the mean radius of the Earth (R E = 6378 Km), we can represent the atmosphere by a stratified 
rectangular slab with the interrogated surface centered on the bottom free. The atmospheric 
extinction coefficient, p (z), is assumed to be a function of the vertical height only and no 
horizontal gradients are assumed. The atmospheric transmission at nadir is then given by 


T^TX/^oc^exp 


-j dz'n{z') 


(Al) 
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where h s is the height of the surface above sea level, and we assume the altimeter is located well 
above the attenuating atmosphere 


The solar irradiance can be viewed as a plane wave incident on the top plane of the atmosphere at 
a solar zenith angle, 0.. If we ignore the slight bending of the solar rays as they propagate ^through 
the vertically stratified atmosphere, the transmitted irradiance to a differential volume at an 
arbitrary height above sea level, z, is given by 


N° exp 


- sec 6 s J dz ' 1 /i(z') = Nl[T(z, qo)] 


sec#. 


(A 2) 


where is the extraterrestrial solar irradiance impinging on the Earth’s atmosphere which, at 
532 nm, has a value N% n = 0.2 watts/m 2 -A°. 

The altimeter receiver has a solid angle field of view (FOV) given by 

(A3) 

where 0 T is the half angle of the far field cone viewed by the receiver. The scattered solar power 
collected by a receiver with collecting aperture A, and spectral bandwidth AX is given by the 
following integral 




N° x {M)A r n r 

An 


( R-zf An 

1 


(A4) 


where (R-z) is the range to the differential scattering volume P(z) is the vohimetnc scattenng 
coefficient at a height z above sea level, and we have assumed isotropic scattermg^Fo 
wavelengths far from absorption lines, the atmospheric extinction is due entirely to scattenng so 
we can set p(z) = p(z). Using (Al) and (A4), the rate at which atmosphenc solar scatter is 
converted into photoelectrons by the receiver is 

. = njh 

K 

oo r r 

J dzju(z) exp -(1 + sec#, ) J dz'/j.{z') 


= 


hv An 

_ %Vr N°{AX)A r Cl 


hv 


An 


(A5) 


where n, is the detector quantum efficiency, r,, is the receiver 

hv is the photon energy at the laser wavelength. In order to evaluate the final integral in (A5), we 
define the new variable 
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JU 

£ = (1 + sec 0,)j dz' p(z') 


(A6a) 


which further implies 

dg = -( 1 + sec 0 s )ju(z)dz 

Substituting (A6) into (A5) and evaluating the trivial integral yields 


(A6b) 


- 


l q *lr 

N‘ t (AA)A,n, 

f 1 Trl+sec#, 

1 " _ 

hv 


[ 1 + sec 0 S 


(A7) 


Note that as T 0 -»1 (no atmospheric scattering), the background count rate correctly goes to zero 
for all values of 0 S . 

When the solar zenith angle equals zero, the background photoelectron count rate due to 
atmospheric scatter has a maximum given by 


and goes to zero as 0 s -Hr/2. 


(A8) 


A.2 SOLAR SCATTER OFF THE SURFACE 

The mean background count rate due to solar scatter off a Lambertian surface is given by the 
equation 

n a = UsOl^L N°{AA)n r pT^ sece ' cos yr (A9) 

7ihv 

where vp is the subtended angle between two vectors - the surface normal and a vector originating 
from the center of the beam pattern on the ground and directed at the Sun as in Figure 3 The 
angle vy can be expressed in terms of the local surface slope, o, and the solar zenith angle, 0 S , via 
the equation 

cosy' = coscrcos#, + sin a sin 0 S cosy) (AlOa) 

where <p is the angle between the plane containing the altimeter line-of-sight and the Sun (x-z 
plane) and the plane containing the altimeter line-of-sight and the surface normal. Since 
2 n > (p > 0 , we have the condition 

cos (a~0 s ) > cosy/ > cos( cr + 0 s ) (A 10b) 
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and cos ^/has the maximum value of 1 when tp-0 and a =0*. Thus, from (A9), the following 
inequality holds 


/?„ 


< 


^^-N%/u)a r pT; 

7U1 V 


(All) 


Note from (A9) that the surface is in shadow (i.e. no scattered surface radiation) when 3rc/2 > vj/ 
> n/2. 


APPENDIX B: EFFECTS OF RECEIVER DEAD TIMES ON SIGNAL CELL 
CONTRAST 

B.l RELEVANT TIME CONSTANTS 

There are several important time constants which can affect the overall performance of the 
microaltimeter receiver. These include: 

t s , the overall spreading of a single laser pulse by the surface under investigation 

the integration time of the photodetector and front aid receiver electronics which 
determines whether or not threshold is achieved in the timing discriminator 

Ta, the “dead time” of the detector or range receiver in which no pulse timing can be 
performed following a triggering signal or noise count above threshold 

t g , the range gate width 

-c b , the duration of a range bin in the correlation range receiver 
Tf, the duration of a flame in the correlation range receiver 
B.2 GENERAL THEORY 

We define the generating function of the Poisson process in the usual way, i.e. 

g(t) = A + s(t) C 01 ) 

where X is a noise background (assumed to be constant over the range gate width) and s(t) is the 
signal waveform which is distributed over a time interval, t s > which varies with surface slope and 
roughness. The mean number of counts occurring in an infinitesimal time interval dt is given by 

dn{t)-g{t)dt 

which is also equal to the probability of a photon being detected in the infinitesimal interval dt by 
a receiver with a detection threshold of one photoelectron . 
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To obtain the true mean photon count within the differential time interval at time t, we must 
multiply (B.2) by the probability that no photon detection took place within the previous time 

interval t-rj <t' < t, i.e. 




J dr g (f). 

t-U 


(B.3) 


Thus, the mean count in the mth range bin, is given by the integral 


mz b 

»m= \ S( 0 ex P 




- j dfg(r) 




dt 


(B.4) 


If signal is not present in the range bin, we can set s(t) =0 in (B.4) and the mean number of noise 
counts collected in a range bin is reduced, relative to a detector/receiver with infinitely fast 
recovery, to 


n b ’ = Xz b exp(-Xr d ) = n b exp(-Xr d ) 


(B.5) 


Thus, for Mhz noise rates and receiver dead times on the order of 100 nsec or less, the mean noise 
count per range bin is reduced by less than 10%. 

Substituting (B. 1) into (B 4) yields the following general result for the mean count within the 
signal bin 


n , ' = exp(-2r li ) 


mz b 

X j dt exp 


(m-l)r* 


i 

- J dt's{t') 




+ 


l 2 

Jd/s(/)ex p 


- } dt'sif) 


t-r d 


(B.6) 


where we have assumed that the signal pulse is confined to the interval t l <t<t 2 which is totally 
contained within the mth range bin. 

The cell contrast in this very general case is now given by 




n k 


n h 


rm b 

X j dt exp 




i 

- f dt's(t') 


t~T d 


+ 


*2 

Jrfrs(/)exp 


J dt's{t') 


(B.7) 


where the subscript d stands for dead time. It can be easily seen that (B.5) through (B.7) 
approach the proper values in the limit as T d -»0, ie. m'- Mib and n,'— Mib + n s and that C approaches 
the value given by Eq. (7a) in the main text. 
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B.3 SPECIAL CASE: Tb^d^. 


If the receiver dead time is longer than the surface return waveform but shorter than the range bin 
(r b > r d >r s =trt,), the second integral in (B.7) can be evaluated by a change of variable, i.e. 

y= = (B8a) 

t-r d 


so that 


*2 

Jr//s(/)exp 


i 

-J dt's(t') 


- = \-e~ n ‘ 


(B.8b) 


To evaluate the first integral in (B.7), we first break the range bin up into four time segments (1- 
before the signal pulse; 2-during the pulse; 3-immediately after the pulse while the detector or 
receiver is recovering; and 4- following detector or receiver recovery). We obtain 
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Using integration by parts [Selby, 1969], we can write for the first integral in (B.9) 
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where we have used the fact that the signal return is a sharp spike within the range bin to pull the 
factor of t out of the integral in (B.10), leaving us with the integral in (B.8b). 

We can perform a similar analysis for the second integral in (B.9) to obtain 
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Substituting (B. 10) and (B 1 1) into (B.9) and inserting the result along with (B8.b) into (B.7) 
gives the following simple expression for the signal cell contrast in the presence of a detector or 
receiver dead time 


c„ = l + (l- 




(B. 12) 
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FIGURE CAPTIONS: 


Figure 1: Normalized surface return rate (fR/fm«x) a function of mean signal strength, n s , and 
three detector thresholds, n, = 1,2, and 3 photoelectrons (pe). Conventional laser altimeters 
typically operate with high signal strengths and high noise thresholds (lower right of the 
asymptotic curve) and therefore make inefficient use of the available laser photons whereas 
photon-counting microaltimeters operate on the most efficient n ( =l curve at low mean signal 

counts (n s <1). 

Figure 2: Principles behind a correlation range receiver. The range gate is divided into equal 
duration range bins, and several consecutive laser fires are combined to form a frame. The 2-D 
area bounded by the range bin and frame dividing lines is a cell. Multiple frames form a 
superframe. Photon counts are accumulated within each cell and, if the total count K exceeds the 
frame threshold K opt , the cell is identified as signal; otherwise it is tentatively identified as noise. 

A valid trajectory is one in which the signal cell is not displaced by more than one range bin in 
consecutive frames. Applying this criteria in algorithms which look forward and backward in 
time and N of M techniques can help recover lost or missing data in near real time. 

Figure 3: Definitions of angles used in the solar background noise analysis. 

Figure 4: Plots of (a) the Differential Cell Count normalized to M, (b) the probability of correctly 
identifying the signal cell in a frame, and (c) the mean number of cells incorrectly identified as 
signal within a frame as a function of the mean signal counts in a frame, N„ for Nb, n ~~ 500 and 
four contrast values, C = 2 (solid), 3 (dot), 5 (dash), and 10 (dot-dash). Note that different 
combinations of C and A, give essentially equivalent performance, but systems with low contrast 
require more signal counts per frame (longer frame times) to achieve the same level of 
discrimination. 

Figure 5: Block diagram of the developmental NASA Airborne Multi-kilohertz Microlaser 
Altimeter Instrument. 

Figure 6: Analysis summary of NASA Airborne Multi-kilohertz Microlaser Altimeter Instrument 
operating at a 10 kHz rate from a cruise altitude of 40,000 ft (12 km): (a) Maximum surface 
sampling rate vs power-aperture product; (b) signal cell contrast vs frame interval, (c) optimum 
frame threshold vs frame interval; (d) normalized differential cell count vs frame interval; (e) 
probability of successfully detecting signal cell vs frame interval; and (f) mean number of falsely 
acquired noise cells per frame vs frame interval. The three curves in each figure correspond to 
different surface reflectances, i.e. p =1 (solid), 0.15 (dot), 0.10 (dash). 

Figure 7: Microaltimeter performance over simulated barren terrain containing a 50 m cliff: (a) 
coarse scale showing the full 4 microsecond range gate; (b) fine scale showing points tentatively 
identified as terrain by the signal extraction algorithm over a short segment on the cliff. 

Figure 8: The following Mars microaltimeter parameters are plotted versus the mean signal count 
per laser fire, n, : (a) single shot probability of detection for 1 pe threshold; (b) required laser fire 
rate to achieve 3 kHz sampling rate from worst case surface, (c) receive telescope diameter; (d) 
transmitted laser energy, (e) average transmitter power; and (f) the Power- Aperture Product. The 
four curves in (c) through (e) correspond to different assumed contrast values, i.e. C =2 (solid), 3 
(dot), 5(dash), and 10 (dot-dash). For each contrast, the corresponding value of the mean signal 
count per frame, N s , is given by Table 5 to ensure equivalent performance. In (f), all four curves 
overlap 
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Figure 9: Block diagram of an all-digital modulo timer integrated with a 6x6 micro-APD array. 
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